This Page Is Inserted by IFW Operations 
and is not a part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of 
the original documents submitted by the applicant. 

Defects in the images may include (but are not limited to): 

• BLACK BORDERS 

• TEXT CUT OFF AT TOP, BOTTOM OR SIDES 

• FADED TEXT 

• ILLEGIBLE TEXT 

• SKEWED/SLANTED IMAGES 

• COLORED PHOTOS 

• BLACK OR VERY BLACK AND WHITE DARK PHOTOS 

• GRAY SCALE DOCUMENTS 



IMAGES ARE BEST AVAILABLE COPY. 



As rescanning documents will not correct images, 
please do not report the images to the 
Image Problem Mailbox. 



MEASUREMENT OF AIR CHARACTERISTICS IN THE LOWER ATMOSPHERE 

RELATED APPLICATION 

5 This application is a continuation of International Application Serial No. 

PCT/AU02/01 129 filed 19 August 2002, published under PCT Article 21(2) in English; and 
claiming priority from Australian patent applications PR 7203 filed 23 August 2001 and PR 
7832, filed 21 September 2001, and applicant claims the benefit of Australian patent 
applications PR 7203 filed 23 August 2001 and PR 7832 filed 21 September 2001 . 

10 

TECHNICAL FIELD 

This invention relates to the use of acoustic signals for atmospheric sounding and is 
particularly concerned with sodar techniques for measuring air velocity variation - such as 
horizontal wind speed variation, wind-shear and/or turbulence - in the lower atmosphere. 
15 The invention may, however, be applied to measuring local density variation in the 
atmosphere, such as may be caused by temperature gradients, temperature, thermal 
inversions and variations in moisture content. 

The apparatus and methods of the intention are also applicable to wind profiling in the vicinity 
20 of airports to enhance air safety and/or permit higher density air traffic at airports. 

The atmospheric sounding techniques of the invention belong to a class of technology 
recently dubbed SODAR, or SOund Direction and Ranging. Sodar is to be distinguished from 
sounding techniques using electromagnetic waves, such as RADAR (RAdio Direction and 

25 Ranging), LIDAR (Light Direction and Ranging), AERI (Atmospheric Emittance Radiance 
Interferometry) and the hybrid RASS (Radio Acousitc Sounding Systems) in the atmosphere. 
However, common to all these techniques in their current form is a concern with Doppler 
signals and the use of Fourier transform methods in processing such signals. While SONAR 
(Sound Navigation and Ranging) has not been mentioned because it is employed in liquid 

30 media, some overlap between the exclusively acoustic techniques of sonar and sodar may 
be seen because sonar ranging and imaging methods have been applied in air - as in some 
camera ranging, non-destructive testing and medical imaging systems. 

BACKGROUND OF THE INVENTION 

35 Though exclusively acoustic methods for wind profiling and the like have a long history, 
Coulter & Kallistratova in their 1999 article "The Role Acoustic Sounding in a High- 
Technology Era" [Meteorol. Atmos. Phys. 71,3-19] show that these methods have not lived 
up to their promise. This appears to have been largely due to an inability to achieve an 
adequate signal-to-noise ratio [s/n]. 



Sodars for atmospheric sounding have almost universally employ short 
(millisecond), single-tone high power pulses, multiple receivers and simple timing 
circuits to determine the sequence of echoes at different receivers needed to 
5 deduce the height of various discontinuities in the atmosphere.. US patent 
2,507,121 to Sivian [1950] disclosed a method for measuring the height of 
atmospheric discontinuities that involved sending such a pulse vertically into the 
atmosphere and, after cessation of the transmitted pulse, detecting vertically 
returned echoes using two similar receivers located near the transmitter. In the 

10 embodiment of most interest, the first receiver is shielded against receiving 
echoes but the second is not and the two receivers are connected so that their 
outputs are opposed and the net signal can be displayed on an oscilloscope. In 
the event of a normally returned echo, a pip is displayed because only the second 
receiver detects a signal. However, In the event of local noise such as a gunshot 

15 both receivers detect the same signal and no pip is displayed. 

US patent 3,889,533 to Balser [1 973] disclosed an 'acoustic wind sensoi* In which 
an acoustic transmitter Illuminates a cylindrical column of air by either CW 
(continuous wave) or pulsed signals and remote receivers are pointed at narrow 

20 or broad portions of two or more sides of the side column to detect acoustic 
energy scattered laterally therefrom. The Doppler components of this scattered 
energy are then extracted to determine wind velocity at the various heights. In 
order to observe a portion of the illuminated column, which is - say - 1000 m 
from the ground, the receivers need to be spaced from the transmitter by a 

25 roughly similar distance. Also, significant spacing is needed to sufficiently 

attenuate the direct signal [sometimes called the 'zero Doppler 1 signal] from the 
transmitter. An application of this system to the detection of persistent vortices 
near runways was disclosed in US patent 3,671 ,927 to Proudian and Balser. 

30 US patent 3,675,1 91 to McAllister [1972] disclosed the use of four adjacent arrays 
of acoustic transducers capable of being used as speakers and microphones, the 
arrays being aligned with the cardinal points of the compass and being shielded 
from one another, except at their upper faces. Short acoustic pulses were 



transmitted vertically upwards and the relative timing of the returned echoes at 
each of the four arrays gave the height and bearing of wind layers. [It might be 
noted that the physics of acoustic sounding was well documented In 1969 by 
McAllister and others In "Acoustic Sounding - A New Approach to the Study of 
5 Atmospheric Structure" In Proc. IEEE Vol. 57, 579-587.] US patent 4,558,594 to 
Balser disclosed the use of an acoustic phased array capable of directing 
successive pulses in different directions, the echoes from one pulse being 
detected by the array before the next is transmitted. US patent 5,521,883 to Fage 
et al uses a similar phased array to send pulses of different frequencies In 
10 different directions and then listen for all echoes simultaneously, thereby 

decreasing the cycle time. The typical angle of elevation for pulse transmission in 
the latter systems was between 20 and 30 degrees. The relatively low elevation 
angle Is to enhance Doppler components in the returned echoes due to horizontal 
(rather than vertical) wind speed in the direction of Interrogation. 

15 

In recent years, radar DSP (digital signal processing) techniques have been 
applied to the sodar to achieve improved s/n. In particular, pulse-compression 
techniques have been used, in which the echoes from a phase or frequency 
coded pulse are processed with matched filters using Fourier transforms to give 

20 the range resolution normally associated with a shorter pulse with a much higher 
peak power. Such coded pulses are said to have 'pulse-compression 1 waveforms 
or to be 'pulse coded'. For simplicity, pulses of this type will be called 'chirps 9 . In 
an article entitled: "Use of Coded Waveforms for SODAR Systems" [Meteorol. 
Atomos. Phys. 71 , 15-23 (1999)], S G Bradley recently reviewed, with simulations, 

25 the use of radar pulse compression techniques to improve amplitude 

discrimination In sodar. Examples of the use of pulse compression techniques In 
radar can be found in US patents 6,208,285 to Burkhardt, 6,087, 981 to Normat 
et al, and 6,040,898 to Mroski et al. Despite the application of such sophisticated 
techniques to sodar, a review by Crescent! entitled, The Degradation of Doppler 

30 Sodar Performance Due to Noise" [Crescentl, G. H., 1998, Atmospheric 

Environment, 32, 1499-1509], found that severe problems remain even at modest 
ranges of 1500 m. 
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OUTLINE OF THE INVENTION 

From one aspect, the Invention comprises methods and apparatus for acoustic 
sounding in air in which echoes from a transmitted chirp (along with extraneous 
acoustic inputs) are detected during transmission of the chirp. In other words, 
5 there fe 'listening while sending 1 . 

This avoids the need to limit pulse length to secure near-range capability, which Is 
essential in known pulsed sodars that employ the transmit then listen' strategy. 
For example, If the pulse length of a conventional sodar is one second, the first 
10 1 70 m of range will be lost because the receiver will be turned off for the first 
second; a 10s pulse will lose the first 1700 m of range. Typically, therefore, 
pulsed sodars of the art employ pulses of a few tens of milliseconds. By contrast, 
our chirps are of at least 300 ms duration and, preferably longer than 10 s; 
indeed, we have used chirps of up to 50 s, the duration only being limited by our 
is current signal processing capacity. Preferably, the duration of the chirp is at least 
5% of the listening time; that is, there is at least 5% overlap between chlip 
transmission and echo reception, but It will be appreciated that listening time 
depends on the distance range covered. For ranges lip to a few km, we prefer 
chirp lengths well over 50% of receive time. As a convenient guide, we listen for 
20 about 6 s longer than the chirp for each km of range. Thus, in a system with a 1 
km range, the chirp/pulse duration might be 15 s and the listening time 21 s; for a 
2 km range, we might use a 31 s chirp and listen for 43 s. Generally, we start 
listening at the commencement of the chirp transmission to obtain data from 
ground level up. For some applications however we may not want the ground 
25 level data and choose to start listening some time after the end of the chirp 
transmission. 

The longer the chirp, the higher its energy for a given transmitter power and the 
better the echoes can be discriminated using appropriate matched filter and/or 
30 Fourier techniques. It is thus much easier to detect faint echoes behind the direct 
signal from the transmitter with long, low power chirps than with conventional 
short high-power pulses. The danger of receiver overload is also mitigated by the 
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use of modern microphones that have a large dynamic range. Of course, acoustic 
shielding of the receivers) from the direct transmitter signal is sensible. 

Indeed, the improvement now possible with the use of long chirps and matched 
filter techniques is such that, from another aspect, the invention comprises 
methods and apparatus for acoustic sounding in air in which the echoes from a 
chirp of greater than 300 ms are detected and processed using matched filter and 
Fourier techniques. Either the 'send then listen' or the listen while sending' 
strategy may be used. As already noted, chirps with duration in the order of 
seconds are preferred; with chirp durations of tens of seconds being favored in 
may situations. 

From another aspect, the invention comprises methods and apparatus in which 
multiple receivers are located near a common transmitter so that each will receive 
echoes from each transmitted chirp. Preferably, the receivers are located close 
enough to share a common acoustic and system noise environment and, 
preferably, they are arranged so as to receive the same direct signal (In both 
frequency spectrum and amplitude). This allows received signal components (eg, 
direct signal, ground clutter and noise) that are common to more than one 
receiver to be to be efficiently removed by differencing the signals from two or 
more receiver locations. Of course, by 'multi-receiver" we mean to include the 
situation where a single receiver is moved to multiple receiver locations and 
where a separate chirp is transmitted for each receiver location. 

While the removal of the common unwanted direct signal, as well as common 
noise components and ground clutter, Is highly desirable, it is very difficult to be 
done directly on the received signals for a chirp that lasts tens of seconds. 
According to another aspect of the Invention, we employ multiple acoustic 
receivers with a single transmitter and process the received acoustic signals in 
the (Fourier) frequency domain using matched filter techniques to generate a 
cumulative phase output for each receiver signal and then manipulate these 
outputs to achieve the appropriate measurement. Subtraction or differencing of 
the cumulative phase signals eliminates the direct signal, common noise, 



common ground clutter and the common slgnaJs due to variation in vertical wind 
speed, the residual differential cumulative phase then represents the variation of 
wind speed with range distance. This overcomes a major problem with 
conventional sodars, which cannot discriminate between returned Doppler signals 
due to vertical wind speed without a direct measurement and those due to 
horizontal wind speed. 

However, comparison or differencing of two or more cumulative phase signals 
requires a common starting or reference point In the signals. This Js conveniently 
the start of chirp transmission, which can be determined by the start of the 
received direct chirp or by an electronic signal from the transmitter. However, 
many other methods of synchronizing the receiver signals are possible. Thus, 
while ft is desirable that the receivers are located in a common acoustic 
environment In the vicinity of the transmitter, It is not essential that they by 
equidistant from the transmitter In order to ensure that the direct chirp arrives at 
each receiver at the same time. 

By transmitting two differently coded chirps (at the same time, using two 
transmitters or one after the other using one transmitter) the cumulative phase 
outputs can be manipulated to remove all common signals, and components due 
to cross-range wind, to allow generation of a further output that is Indicative of 
variation of the speed of sound with range and, thus, variation of temperature with 
range. Preferably, the two chirps are Identical positive and negative linear phase 
chirps (eg, the positive one rising from 800 to 1600 Hz and the negative one 
descending from 1 600 to 800 Hz at the same phase rate. 

Thus, the last-mentioned aspect of the invention provides a further large 
improvement in s/n, allowing much Improved echo discrimination with respect to 
the art, despite listening while sending. Also, simultaneous echo reception and 
processing by multiple receivers greatly Improves cycle time. 

A convenient arrangement of receivers In a system for vertical atmospheric 
sounding is to locate one receh/er at each cardinal compass point around the 
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transmitter and to slightly Incline opposed receivers toward or away from one 
another. Thus, the phase components common to the N-S receiver signals are 
removed by phase differencing to leave that associated with variation of the net 
N-S wind over range distance. Systems of this type are suitable for vertical 
5 sounding in noisy environments such as airports, power stations and urban areas. 

It will be appreciated that the invention is not limited to the use of four receivers, 
or to vertical sounding systems or to the symmetrical placement of receivers 
around a transmitter. The receivers may be arranged in a line, for example across 

10 a runway glide path with one or more transmitters to detect persistent vortices 
caused by the passage of large aircraft. The high-speed localized winds which 
can make up such vortices are difficult to quantify because the high DoppJer shifts 
of echoes generated have considerable ambiguity. In this situation, another 
aspect of the invention involves the repeated analysis of recorded echo signals 

15 for each range point using a matched filter that is fed with a succession of 

different reference chirps so -that, for each range point, a reference chirp is found 
that generates a zero phase gradient output (In the Fourier domain). That 
reference chirp is then indicative of the wind speed at that range point. It is also 
useful here (as well as in other applications envisaged by the Invention) to take 

20 'readings* In the absence of vortices to record ambtent wind and noise conditions 
and to subtract the associated phase signals from those generated when there is 
a vortex present 

Whilst listening during sending is not essential for the implementation of the last 
25 described aspects of the invention, it is certainly desirable because It enables the 
use of long chirps, better echo discrimination and the effective elimination of 
range dead-zones. 

The transmitted acoustic chirp can be generated by feeding a commercially 
30 available loudspeaker (transmitting acoustic transducer) with an electrical Input 
signal from the sound card of a computer (for example), while the echoes can be 
detected using commercially available microphones (receiving acoustic 
transducers). The loudspeaker and microphone(s) can be mounted with separate 



concentrating reflectors (plates, horns, dishes or the like) or they may be mounted 
with a common reflector. For example, four microphones can be arranged in 
quadrature around a single loudspeaker above a single reflector dish so that the 
transmission axis Is substantially axial with respect to the dish. Since the 
microphones are then offset with respect to the axis of the dish, the receiving 
axes of opposed pairs will be oppositely inclined toward the transmission axis; 
that ls r each receiver will be most sensitive to echoes coming from a direction 
opposite to its location oh the dish with respect to the transmission axis. 

An arrangement where a loudspeaker and multiple microphones are mounted on 
a common structure allows the transmission axis to be conveniently aimed or set 
as desired by moving the structure. A system with a transmission axis of low 
elevation can be used, for example, to detect or characterize vortices caused by 
large aircraft landing or taking off at an airport It Is even possible to use airborne 
systems of this type to warn pilots of clear air turbulence (CAT) that is difficult to 
detect using radar. For example, a compact transmitter and receiver system could 
be mounted in the nosecone of a large aircraft. Alternatively, only the transmitter 
need be mounted in the nosecone since the receivers can be mounted in a row 
along the leading edge of the wings. 

As already noted, the receiving axis of a receiver may be inclined with respect to 
the transmission axis and that, where multiple receivers and signal differencing 
are used, it is desirable that the axes of opposed receivers are equally inclined. 
The optimum angle of inclination will depend upon the aperture of the receivers, 
the range of the system and the desirability of beating the receivers in a common 
acoustic environment. An angle of 20 degrees in a system with a 3 km range is 
likely to place the receivers too far from one another to have a common acoustic 
environment, but this may not be so In a system with a range of only 250 m. 
Angles of Inclination of between about 2 and 10 degrees have been found 
suitable, with angles between 4 and 7 degrees preferred. It will be seen that the 
point of Intersection of a receiver axis with the transmission axis is not intended to 
be the nominal range of the sodar system. Indeed, highly satisfactory results have 
been obtained where medium aperture microphones are located about 1 m from 
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the loudspeaker in a common dish with their receiving axes angled at about 4 
degrees to the transmission axis. In effect, the receivers of an opposed pair are 
looking for wind-Induced Doppler signals from large Illuminated areas on opposite 
sides of the transmission axis but in the same plane as the receiving axes and the 
transmission axis. 

As already noted, it is desirable (but not necessary) to space multiple receivers 
equidistant from and near to a common transmitter so that each will be subject to 
the same ambient noise (as well as other common components). Generally, the 
louder and less uniform the noise environment, the nearer the receivers need to 
be to one another to ensure that each is subjected to the same environmental 
noise, as far as practicable. We have found that, in a noisy environment, the 
distance between a receiver and the transmitter should be of the order of metres. 
In a quiet environment, it can be of the order of 1 0 m. 

In general, the transmitted chirp should have a tonal range (acoustical bandwidth) 
suited to the object being sounded. We have found that wind-shear below 3000 m 
is best sounded at the lower end of the audible range; for example, 500 - 5000 
Hz, more preferably between 800 Hz and 3 kHz and most preferably between 1.0 
kHz and 2.5 kHz. 

While the tones In a chirp can be modulated in frequency and/or phase In many 
ways in conformity with pulse compression techniques, we have found It 
convenient to use a linear chirp In which the frequency increases or decreases 
smoothly and linearly from start to finish of the chirp. Ideally, such a chirp has a 
uniform rate of phase-shift. The use of positive and negative linear chirps Is of 
particular value in the reduction of unwanted phase components in techniques for 
measuring air temperature disclosed herein. Linear chirps are also easily 
generated and their echoes convenient to process using available DSP and 
Fourier techniques implemented using personal computers. • 

Although (as already noted) long duration chirps offer the potential of high system 
processing gains (lower s/n), long chirps also result in significant computational 
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demands when using the high signal sampling rates and the Fourier techniques 
needed to achieve such gains. We have found that current readily available FFT 
algorithms, DSP chips and PCs set a practical limit on chirp duration of about 40 
- 50 s at sampling rates of about 96 k Hz. This typically represents some 1400 
samples per m, given a range of 3000 m. Indeed, the computational demands are 
such that we prefer to dedicate one PC to each receiver of a multi-receiver 
system so that echo analysis for all receiver signals can proceed in parallel to the 
point where signal differencing takes place. In the future, developments in chips, 
FFT/matched filter techniques and PCs may allow longer chirps to be processed 
using a single PC - or, much faster updating times using the pulse lengths 
presently achievable. 

DESCRIPTION OF EXAMPLES 

Having portrayed the nature of the present invention, particular examples will now 
be described with reference to the accompanying drawings. However, those 
skilled in the art will appreciate that many variations and modifications can be 
made to the chosen examples while conforming to the scope of the invention as 
outlined above. 

Brief Description of the Drawings 
In the accompanying drawings: 

Figure 1 is a series of diagrammatic plan views showing selected 
arrangements of transmitters and receivers, the transmitters (loudspeakers) being 
shown as small shaded circles and the receivers (microphones) being shown as 
small unshaded circles. 

Figure 2 is a series of diagrammatic elevations showing co-located and 
separately located transmitter and receiver arrangements. 

Figure 3 is a diagrammatic sectional elevation showing the arrangement 
of the transmitter and receivers of the first system example. 

Figure 4 is a schematic plan of the system of Figure 3 showing the 
genera] manner in which signals to the transmitter are generated and signals from 
the receivers are processed. 
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Figure 5 is a block diagram illustrating a system for extracting east-west 
phase information from the echoes received by the east and west microphones of 
the chosen system example. 

Figures 5A, 5B and 5C are block diagrams Illustrating respective parts of 
» a circuit and process by which wind-speed, wind-bearing and wind-shear 
information with respect to height can be generated with the system of the first 
example. 

Figure 6 is a histogram depicting a typical digitized acoustic signal 
detected by a receiver of the system of the first example. 

Figures 7 A and 7B are graphs of typical signals before and after the tow- 
pass filter of the matched filter of the system of the first example. 

Figure 8 Is a series of graphs depicting the cumulative phase information 
derived from all receivers of the system of the first example. 

Figure 9 Includes bar charts and graphs showing the atmospheric wind 
characteristics output from the system of the first example. 

Rgure 10 is a block diagram illustrating portion of the process by which 
temperature with respect to height is derived 

Figure 11 is a graphical representation of how phase signals of positive 
and negative chirps are differenced to yield an indication of temperature with 
respect to range distance. 

Figure 1 2 is a plot of temperature with respect to range distance 
generated by the use of the disclosed system. 

Figure 13 is a block diagram Illustrating portion of the process by which 
vertical wind speed with respect to height is derived. 

Rgure 1 4 Is a schematic diagram illustrating the detection of wake 
vortices left by a large airplane near a runway. 

Rgure 15 is a graph showing the notional variation of excess phase with 
respect to range distance, this relatlonshp being used to compute wind velocity 
and height In awake vortex. 
Description of the Preferred Brfcodiment 

As illustrated by the simple plan-view diagrams (a) to (I) of Figure 1 f there are 
many possible ways in which the transmitter and receivers of sodar systems 
envisaged by this invention can be configured, and many others are possible. 
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Diagram (a) shows a convenient and economical configuration in which four 
microphones 10 are spaced within and around a common parabolic reflector dish 
1 2 and a single loudspeaker 1 4 is located at the central focus of dish 1 2. In this 
way, the amplitude lobe of the transmitted pulse is vertical but the amplitude lobe 
s each received echo is angled slightly to towards the axis of dish 1 2 and 
loudspeaker 14. 

Diagram (b) of Figure 1 shows three microphones 16 evenly spaced in a common 
dish 1 8 that also has a centrally located loudspeaker 20. Diagram (c) shows a 
io dish 22 with a central loudspeaker 24 and only one offset microphone 26, dish 22 
being mounted so that It can be rotated to successively put microphone 24 in 
different positions [eg, those, illustrated In (a) or (b)]. Diagram (d) shows four 
microphones 28 mounted in a common receiving dish 30 that is separate and 
spaced from the associated loudspeaker 32, which is mounted In its own 
15 transmission dish 34. Optionally, receiving dish 30 may have more or less than 
four microphones located therein. Diagram (e) illustrates a configuration In which 
each of four microphones 36N, 36S, 36E and 36W has Its own receiving dish 
38N, 38S, 38E and 38W (respectively) and a single loudspeaker 40 has its 
separate dish 42. Diagram (f) is similar to <e) except only three microphones 44, 
20 46 and 48 and their respective dishes 50, 52 arxl 54 are deployed around a single 
microphone 56 and its dish 58. Finally diagram (g) shows a configuration in which 
a single microphone 60 and its dish 62 are mounted so as to be rotatable around 
a single loudspeaker 64 and Its dish 66, so as to be able to simulate 
configurations such as those of (e) and (f). 

25 

Figure 1(h) is a plan view of a linear array of four pairs of receivers 70a and 70b, 
71a and 71b, 72a and 72b. and, 73a and 73b arranged In a row with one receiver 
of each pair positioned on either side of a single transmitter 74. Transmitter 74 
generates a narrow spherically propagated beam orthogonal to the line of 
30 receivers so that the signals received by each receiver pair can be processed to 
remove the direct signal and common noise. In the linear array of Figure 1(1), 
which is also shown in plan view, A row of transmitters 76 is arranged parallel to a 
row of receivers 78, a common signal being fed to all transmitters so as to 
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generate a linearly propagated sound wave. The arrangements of Figures 1h) 
and 1 (Q are suited for arrangement across the glide path of an airport to detect 
persistent vortices. [A system of this nature will be described in more detail 
below.] 

The diagrams (A), (C) and (E) of Figure 2 are diagrammatic elevations of 
configurations (a), (c) and <e) respectively of Figure 1 . 

The system 100 of the first example, illustrated In Figure 3, approximates the 
arrangement of Figure 1 (a) and Figure 2(A) and Is adapted for vertical or inclined 
atmospheric sounding where both cross-range (horizontal, in this case) and 
along-range (vertical, In this case) wind velocities are required. For convenience, 
however, it will be assumed that chirps from the transmitter are directed vertically 
upwards to illuminate an inverted cone of air indicated at 101 by broken lines. 
Also for convenience, the cardinal points of the compass will be referred to as N, 
S, E & W, as well as north, south, east and west where thought necessary. 

System 100 includes a large main dish 102 and a small secondary dish 104 
mounted directly above the main dish. A transmitter/receiver module 106 is 
supported centrally above large dish 1 02 (by struts that are not shown) and, in 
turn, supports small dish 104 on the top thereof. Module 106 comprises a sound- 
adsorbent molding 108, Into the bottom of which a central loudspeaker 110 and 
four peripheral microphones are fitted. 

The microphones are arranged in quadrature and aligned with the cardinal points 
of the compass so that, in the sectional diagram of Figure 3, the W microphone is 
shown at 112 (on the easVright side of loudspeaker 1 10) and the E microphone is 
shown at 114 (on the west/left side of loudspeaker 110). This apparent reversal of 
naming the E and W microphones is convenient because the microphone on the 
west side of the loudspeaker Is positioned to be most sensitive to echoes coming 
from the east, after reflection and focusing by dish 104, and vise versa. The axis 
of each microphone is angled to the vertical at between about 3 and 10 degrees. 
The loudspeaker 1 10 and microphones 1 12 and 1 14 are located near the focus of 
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large dish 102. A fifth directional microphone 1 16 Is located at the focus of small 
dish 104. 

In the diagram of Figure 3, a single horizontal reflective atmospheric discontinuity 
5 (such as the nocturnal boundary layer or other thermal inversion layer - TIL), 1 20 
Is shown. Since loudspeaker 1 1 0 is pointed vertically downward, ft will generate a 
downwardly directed vertical beam that will be reflected vertically upward by large 
dish 102 along a central system axis 122. Some echoes will be returned down 
axis 1 22 to microphone 1 1 6 on small dish 1 04. However, the beam of 

10 interrogating pulses will be conical and will Illuminate a significant area of the TIL 
120 around axis 122 and echoes will be returned from an area to the west of axis 
122 along the axis 124 of west microphone 1 1 2 and be most strongiy detected by 
that microphone (In comparison with the signals detected by the other 
microphones). Similarly, echoes from TIL 120 to the east of axis 122 will travel 

is along the axis 126 and be most strongly received by east microphone 1 14. 

While echoes returned from TIL 120 and detected by W microphone 1 12 may be 
centered about path 124, microphone 1 12 will pickup echoes from a large area of 
TIL 120 in the vicinity of axis 122. Thus echoes from a source near axis 122 are 

20 likely to be picked up by all microphones. Doppler (phase) components common 
to echoes detected by all microphones at much the same time are therefore 
Indicative of the vertical velocity of TIL 1 20 in the vicinity of axis 122, and it can be 
expected that these common components will be most prominent in the echoes 
detected by central microphone 1 1 6. If the Doppler components of echoes 

25 received by W and E microphones 1 1 2 and 1 1 4 are subtracted, the common 
Doppler (phase) components indicative of the vertical velocity component will be 
removed and It can be assumed that the remaining Doppler (phase) components 
are due to net wind speed in the east/west direction. Similar subtraction of the 
Doppler components of echoes received by the N and S microphones will yield. 
30 the net wind speed in the north/south direction. 



In practice, of course, there will be many atmospheric discontinuities at many 
altitudes within range that generate echoes and that the time of return of such 
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echoes will be Indicative of range or attitude and the amplitude of the echoes will 
be indicative of the magnitude of the respective discontinuities. 

The generation of chirps for transmission and the processing of received echoes 
may be implemented in many ways, whether in software or hardware. The mode 
of implementation will be influenced by the desired chirp length, listening 
time/intended range, sampling rates, and up-date frequency, since these factors 
largely determine computation demand. To provide desired high processing gain, 
chirp durations greater than 0.3 s are considered essential, with durations of tens 
of seconds desirable. Prior art pulsed systems using the transmrt-then-llsten 
strategy typically have transmit times In the order of tens of milliseconds and a net 
listening time of about 6 s for a range of 1000 m. By contrast, in the present 
example, for a range of 1000 m and net listening time of 6s, the selected chirp 
duration is 37s and the total listening time is 43 s. For the same total transmitted 
energy the chirp of the present example can have a thousand-fold lower peak 
power than a 37 ms chirp typical of the prior art The total listening time of the 
present example is more than seven times that typical of the art for a 1000 m 
range, providing much greater opportunity for processing gain. By using chirps 
(pulse-compression waveforms) and matched filter processing, processing gains 
are further multiplied many-fold. 

The computation demands of the system of the example are, however, 
substantial and, In this case, were thought to Justify dedicating a PC to process 
the signals from each receiver and using another PC as a controller. Figure 4 
shows this arrangement In which the N and S receiver microphones are shown at 
130 and 132 respectively, the dedicated N, E, W & S and vertical (V) calculating 
PCs are shown at 134, 136, 138, 140, 142 and 143 respectively and the controller 
PC is shown at 144. Controller PC 144 generates the chirp for transmission by 
transmitter/loudspeaker 1 10 and the reference chirp for use In matched filtering 
by the calculating PCs. It also collects the results of the computations of PCs 134 
- 1 42 for integration, display and reporting. 
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In this example, the chirp has a phase/frequency that increases linearly over the 
37 s from 800 Hz to 1600 Hz (the chirp could just as easily decrease linearly from 
1600Hz to 800Hz) and is emitted at an acoustic power of a few hundred milli- 
watts that remains constant for the duration of the chirp. This type of pulse 
s compression waveform has a small bandwidth (about 800 Hz) and Is simple to 
generate accurately using a PC sound card and conventional loudspeaker driver 
circuits for powers up to many Watts. It is also one for which a 'matched fitter' or 
correlator can be readily designed and used to extract echoes of the china from 
received signals having high noise levels and to effectively compress the energy 
w of each echo Into a period of time that is much shorter than that of the transmitted 
chirp. The use of pulse-compression waveforms and matched filters thus yield 
very high processing gains. 

Measurem ent of Cross-Ranae Wind Velocity 

15 System 100 of the chosen example Is well suited for measuring various 

characteristics of horizontal (cross-range) wind. The manner in which the net 
phase difference between the east and west echoes is effected to derive the E-W 
wind speed with range distance will now be described with reference to Rgures 
5A, 5B and 5C. H should be appreciated that the system of Figure 5A is 

20 duplicated for the extraction of the net phase difference between the north and 
south echoes in exactly the same manner. This Is Indicated by WN" and "E/S" in 
Figure 5A. However, the following description will refer principally to the 
processing of the E & W echoes. 

25 The processes of Figure 5A are performed in dedicated E & W PCs 138 and 140 
and control PC 144 to (i) extract the cumulative phase and amplitude Information 
from the noisy received signal using matched filtering and (ii) difference the E - W 
cumulative phase information to eliminate the zero Doppler signal (and other 
common noise components) to yield a Doppler phase signals indicative of net 

30 east - west wind speed. The cumulative phase information and the differenced E 
- W and N - S phase information are illustrated graphically In Figure 8. The E - 
W and N - S signals are then used as inputs for the circuit and process of Figure 
5B, from which overall wind velocity and wind shear information Is derived. PCs 
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138 and 140 again do the processing, the results being collated and displayed on 
control PC 144 and illustrated graphically in Figure 9. Figure 5C illustrates in more 
detail the phase and amplitude extraction circuit, which includes a matched filter 
or correlator, shown in Figure 5A. The phase extractor Is substantially identical for 
each receiver and Is implemented in the respective receiver PC. The circuit and 
process of Figure 5 are illustrated and described generally because the 
techniques employed are analogous to known signal processing techniques used 
in radar. 

In Figure 5A, loudspeaker 1 10 Is shown pointing upwards, the transmitted chirp is 
indicated by shaded graph 150, east microphone 1 12 Is shown pointing up on the 
right side of the diagram to receive input signals Indicated by arrow 152 and west 
microphone 1 14 is shown pointing up on the left side of the diagram to receive 
input signals indicated by arrow 154. Chirp 150 is generated using a control PC 
output 156 (indicated by associated graph 169) to drive a transmitter circuit 158, 
which in turn, powers matched loudspeaker 1 10 to ensure that the acoustic power 
of transmitted chirp 150 is uniform over its whole bandwidth (about 800 Hz). In 
this example, the power is about 0.5 Watt; not enough to cause annoyance even 
In urban areas. 

Received east and west signals 152 and 154 are respectively processed 
(essentially amplified and band-pass filtered) in receivers 159 and 160, the 
outputs of which are sampled at 96 kfe and digitized by analog-digital (A-D) 
circuits 161 and 162 to generate digital receiver signals 163 and 165. A 
representation of signal 1 63 or 1 65 is shown in Figure 6. 

Outputs 163 and 165 are fed to phase and amplitude extraction circuits 166 and 
167 that employ matched filters to correlate signals 163 and 165 in the Fourier 
domain with a version of the transmitted chirp, extract or derive echo phase and 
amplitude information therefrom. Said version of the transmitted chiip is provided 
via line 168 by control PC 144 to each extractor 166/167 and comprises the 
transmitted chirp shifted down by 800Hz; ie, converted to a 0-800 Hz chirp as 
indicated In nearby graph 169. It may also be desirable to delay signal on line 168 
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with respect to the actual transmitted signal so as to select a desired range band 
of the system. 

Extractor circuits 166 and 167 each have two outputs. Outputs 170 and 171 of 
extractor 166 are respectively indicative of the cumulative phase and amplitude of 
the Input 163 after matched filtering (Fourier processing), the east cumulative 
phase being shown in Figure 8(l), the south cumulative phase being shown in 
Figure 8(v) and the east and south magnitudes being shown in the barchart of 
Figure 9(l). Outputs 172 and 173 of extractor 167 are respectively Indicative of the 
cumulative phase and amplitude of the input 165 after matched filtering, the west 
cumulative phase being shown in Figure 8(H), the north cumulative phase being 

shown in Figure 8(lv) and the west and north amplitudes being shown In the 
barchart of Figure 9(1). Cumulative phase outputs 170 and 172 contain the 
information from which wind velocity, direction and wind shear can be derived, as 
well as common phase noise and zero-Doppier components due to the direct 
signal and ground clutter. Cumulative phase outputs 170 and 172 are differenced 
in circuit 174 to remove the common components and output on line 176 as the 
net E-W phase or cumulative phase of graph (iii) of Rgure 8. These outputs 
contain the information from which the E-W velocity and wind shear are derived 

As already noted, the circuit and process of Figure 6 is applied in an Identical 
manner to the north and south signals to derive the cumulative north phase 
illustrated graphically by Rgure 8 (Iv), the cumulative south phase illustrated by 
Rgure 8(v) and the net north-south phase variation with height Illustrated by 
Rgure 8(vi). For convenience, east-west phase difference output from circuit will 
be Identified as 176 E-W of and the north-south phase difference will be identified 
as 176 N-S. 

Rgure 5B illustrates how outputs 176 E-W and 176 N-S are processed in the 
chosen example to generate wind speed and bearing (ie, wind velocity) and wind 
shear magnitude with range distance (altitude, in this case). Outputs 176 E-W 
and 1 76 N-S are fed to a combiner circuit 1 80 from which wind bearing or 
direction Is derived, as shown in Rgure 9 (II). This is done by determining whether 
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the E-W signal Is positive or negative and whether the N-S signal is positive or 
negative in order to place the wind direction in the correct quadrant of the 
compass. Thus, a positive E-W phase and a positive N-S phase will indicate that 
the wind direction is In the first quadrant, and so on. The angle In the first 
5 quadrant can be determined by computing the relevant vector to give a more 
precise indication of bearing. 

By employing the well known square-root of the sum of the squares algorithm, the 
magnitude of the phase signal can be determined to output the wind speed and 
10 wind shear magnitudes. This Js done by feeding outputs 176 E-W and 176 N-S to 
respective squaring circuits 182 and 1 84, summing the outputs of these circuits in 
adder 186 and deriving the square root in circuit 188.The resultant output on line 
1 90 is Indicative of the variation of wind speed with respect to altitude. This is 
illustrated by Figure 9 (111). 

By taking the derivative of the signal on line 190 using differentiator circuit 1 92 the 
magnitude of wind shear with altitude is output on line 194. Figure 9 (rV) illustrates 
this. Signal on line 194 can be further processed to apportion the wind shear into 
various magnitude bins using a differencing circuit 196 to generate a series of 
outputs 1 98 - 202 etc, which are illustrated by the graphs of Figure 9(v). 

Referring now to Figure 5C, the basic operation of each extractor will now be 
described. Since the extractors are substantially identical, only extractor 166 for 
the east receiver will be described. Again it will be appreciated that the matched 
filter and the phase and amplitude extraction techniques adopted here are known 
to those skilled in the radar art Preferably, as also mentioned before, each 
extractor is implemented In a separate PC that is dedicated to one receiver. 

Extractor 166 essentially comprises a matched filter 210 that - in the Fourier or 
frequency domain - matches a reference chins (usually a version of the 
transmitted chirp) to the confused time-domain echoes that are included in input 
signal 163. Each digitized sample of E input 163 is converted to complex form. 
The Imaginary part (I) is generated by mixing the input with a digitized 2000 Hz 
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sine signal 300 supplied by control PC 144 using multiplier 302, and the real part 
(Q) Is generated by mixing the Input with a digitized 2000 cosine signal 304 (also 
supplied by PC 144) using multiplier 306. The resultant I and Q signals for every 
sample taken during the listening period are fed to a complex FFT [fast Fourier 
transform] process 308, where all samples are presented and processed as an 
array to generate the frequency domain outputs P and Q\ These outputs are low- 
pass filtered at 310 to remove the upper side band. Figures 7A and 7B show the 
unfiltered and filtered signal components Q' and Q", respectively. The filtered 
signals I" and Q" are then passed to complex multiplier 312 in which they are 
multiplied with the down-converted output signal 168 from the control PC 144. 
The result of this multiplication is then subjected to complex Inverse FFT at 314 to 
generate real and imaginary sample-like outputs I"' and Q'" In the time domain. 

Each 'sample' output I" and Q" from JFFT 314 Is then processed to provide 
corresponding phase and amplitude outputs 170 and 171. Phase output 170 is 
generated by implementing the function Atan2(IYQ") in process 318 to yield a 
succession of phase values between -n and +k on line 318, which is input into an 
unwrap process 320. The unwrap process is known and implementations are 
available in programs such as MatLab. Essentially, process 316 counts the 
number of 2n phase shifts to generate an accumulated phase. For every 
transition from +k to -w (an increasing phase) the phase accumulator is increased 
by 1, and vice versa. This output can then be displayed as a radian count with 
respect to sample number (proxy for time and distance) and displayed graphically 
as in Figure 8. The amplitude output 171 is generated by implementing the 
function [I 2 + Q a ] 1/2 in process 322. The barchart (I) of Figure 9 shows the 
variation of amplitude with height, each amplitude reading being color-coded to 
indicate magnitude. It will be appreciated that a sharp variation in processed 
signal amplitude at a given height is indicative of moisture or temperature change 
and not of wind shear. However, significant temperature differentials In the 
atmosphere will be accompanied by wind change. 
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Measurement of Air Temperature 

Using the system of the above example, the variation of air temperature over 
range distance can be estimated by the use of positive and negative chirps and 
manipulating the cumulative phase outputs generated. This comprises the second 
detailed example of the application of this invention. 

Though not essential, the use of substantially identical positive and negative 
linear chirps is highly desirable in the means for the measurement of temperature. 
Consistent with the above example, a positive chirp that rises in frequency from 
800 to1600 Hz and a negative chirp that falls in frequency from 1600 to 800 Hz 
over a period of 37 s will be assumed. It will also be assumed that the positive 
chirp is transmitted first and that the negative chirp is transmitted immediately 
after the listening time of 43 s has elapsed, there also being a listening time of 43 
s after the negative chirp has been transmitted. In this case, it will be convenient 
to separately digitize the acoustic signals from each of the four receivers (N, S, E 
& W) for the positive chirp and for the negative chirp and to then process each in 
the manner describe above to extract the respective cumulative phase signal. 

The manner in which the cumulative phase signals are manipulated will now be 
described with reference to Figure 10, In which the N, S, E & W cumulative phase 
signals for the positive chirp are shown as inputs 300 and the N, S, E & W 
cumulative phase signals for the negative chirp are shown as inputs 302. Inputs 
300 are added together and divided by four In adder 304 to remove the horizontal 
(cross-range) wind components for the positive chirp, and inputs 302 are added 
together and divided by four In adder 306 to remove the horizontal (cross-range) 
wind components for the negative chirp. The use of four receivers for temperature 
measurement is convenient as the same system can be used for horizontal wind 
measurements as well. The temperature measurement can also be made by 
using a single receiver pointed vertically but such a system could not be used for 
horizontal wind measurements. The outputs of adders 304 and 306 are then 
differenced In process 308 to remove common components due to direct signal, 
vertical wind, ground clutter and noise, leaving temperature related cumulative 
phase difference. The gradient of this difference signal is then derived In process 
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310 and normalized In process 312, the output of this process being a measure of 
the relative change of temperature with altitude (range). To calibrate this, the 
actual temperature near ground level is Input at 314 and a chart - indicated at 
316 - of temperature variation with altitude can be generated. An actual chart 
generated by the means of the second example is appended as Figure 

The physical basis of temperature measurement In the system of this example is 
that the total rate of phase advance with respect to distance of a positive chirp 
passing through a layer of cold air will be slightly less than that for a negative 
chirp, the difference being dependent upon temperature. The total rate of phase 
advance is made up of the sum of a dominant component due to the propagation 
of sound in air (nominally 1 4,500* radians for a distance of 1 km) and a minor 
component due to the internal rate of phase advance within the chirp ( eg, +800 x 
2*743 or 18 radfans/s for a positive chirp and -800 x 2n /43 or 18 radians/s for a 
negative chirp). In the case of a positive chirp, the internal rate is positive and is 
added to the propagation rate; with a negative chip the internal rate is negative 
and is subtracted from the propagation rate. Thus, the rate of Increase of 
cumulative phase with respect to distance is slightly less for the negative chirp 
than for the positive chirp. However, when a cold layer of a fixed distance is 
encountered, the propagation rate of the positive chirp is slowed slightly so that 
the chirp takes longer to travel the distance the internal phase advance is 
proportionately greater than it was when traveling the same distance in warmer 
air. And, since the increased Internal phase advance is positive, it will add slightly 
to the (now slower) propagation^ phase advance. While the cold layer also slows 
the propagation (phase rate) of the negative chirp and the internal (negative) 
phase advance is also increased as a result, the marginal increase Is not as great 
due to the slower phase advance of the negative chirp and therefore results In a 
smaller total cumulative rate of phase change, resulting in a divergence of the 
rate of phase change with respect to time/distance. This is illustrated in the 
diagram of Figure 11. 
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Measurement of Down-Ranoe Wind Velocity 

There are two ways of estimating down-range (in this case vertical) wind in 
accordance with the principles disclosed herein. The first (comprising the third 
example) is to use the central receiver 1 16 in a send-then-listen mode; the 
5 second (comprising the fourth example) is to use the four N, S, E &W receivers in 
a listen-while-sending mode. In both cases, long chirps (greater than 300 ms) are 
employed as taught herein but, in the first, there is some sacrifice of low-attitude 
range and, in the second, there is some sacrifice of accuracy. 

10 in the third example, a short chirp of about 0.5 s Is sent by loudspeaker 110 and, 
immediately after, the incoming signals to microphone 116 are processed by the 
phase and amplitude extraction process described above to generate an 
corresponding cumulative phase and amplitude output, from which vertical wind 
speed variation can be read or deduced. This method will lose the first 85 m of 

15 range and will be subject to errors due to noise and cross-range wind speed, but 
the greater resolution offered by the long chirp will be gained. 

In the fourth example, a chirp of about 5 s can be employed while all receivers 
are listening, the acoustic outputs of each receiver then being processed and the 

20 respective phases and amplitudes extracted as described above. Referring to 
Figure 13, the cumulative phases of the N, S, E & W receivers, indicated at 320, 
are fed to adder 322 where they are added together and divided by four, as in the 
temperature measurement case to remove common elements due to cross-range 
wind. System dependent phase shift 323 is then removed In process 324, tt being 

25 easily removed because it has a constant gradient The output of process 324 Is 
indicative of the variation of down-range wind speed, but is degraded by the 
(relatively short) direct signal and by ground clutter. These undesired signal 
components are manifest in a relatively large phase signal at the origin (zero 
distance and time) Since ft is reasonable to assume that the wind speed at 

30 ground level is zero the initial or residual phase can be subtracted from the 
cumulative phase in process 326 and normalized in process 328, yielding the 
display indicated at 330. Alternatively, the resultant phase can be scaled and 
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calibrated according to a known wind speed at a given altitude obtained by other 
means, such as radio sonde. 

Detection of Vortices Near Runways 
5 In this fourth and final example, an array 400 of receivers and transmitters - such 
as indicated in Figure 1 (h) or (i) - is arranged across the glide path at the end of 
an airport runway. Figure 14 shows the use of an array of the type shown in 
Figure 1 (h) having a single central transmitter 402 and multiple receivers 404 
extending in a row on each side. The entire array may span 150- 200 m. In 

10 Figure 14 transmitter 402 is shown generating a spherically propagated acoustic 
chirp 406 and the wake vortices are shown at 408, arrows 410 indicating the 
direction of rotation of each vortex. The Interaction of chirp 406 with the vortices 
results in the backscatter of echoes 412, which are picked up by receivers 404 
and fed to respective matched filters in a manner similar to that described in the 

is first example. Multiple receivers are used here to provide better horizontal 
resolution and increase overall receiver gain. 

The difficulty in this case Is, however, that the high Doppler echoes returned from 
the vortices creates considerable ambiguity in the results so that the location, size 

20 and speed of the vortices cannot be measured with sufficient accuracy by using 
the system of the first example alone. It is necessary to obtain measurements of 
the ambient conditions prior to the arrival of a large plane and use those 
measurements to adjust and sharpen those taken with a vortex present, after the 
plane has passed. The ambient measurement provides a reference for the 

25 amplitude and phase of the system as well as for the ambient amplitude and 
phase. These results are stored. 

When the echoes from a vortex are being processed, the signal processing 
proceeds for each receiver as described In the first example, except that the 
30 return signal is correlated against multiple different Doppler shifts (positive, zero 
and negative) for the internal multiply of the matched filter using a corresponding 
series of reference chirps generated by the control PC. Each phase calculated In 
respect of the vortex condition for each receiver is then differenced with the 
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corresponding calculated phase for the ambient condition for each respective 
receiver, removing system and noise phase shifts. This leaves a residual or 
'excess* phase shift that can be plotted with respect to distance (time) along with 
the corresponding Doppler shifts, as shown in Figure 15. 

To then estimate Doppler vs. distance (height), it is necessary to search through 
each excess phase record for each Doppler shifted result to find the distance 
ranges for which the excess phase is zero. By combining the ranges for which a 
zero excess phase is found a graph of Doppler shift vs. distance can be made. 
From this it is an easy matter to estimate wind speed vs. distance. It is to be 
noted that the excess phase alone cannot be used to without the additional 
Doppler processing because the range ambiguity effect spreads the excess 
phase result out to the extent that the location and size of the vortex cannot be 
sufficiently located. 

Bradley (cited above) discusses a similar estimation process for the optimization 
of amplitude and reference should be made to that paper. However, using phase 
is more accurate and stable than amplitude, allowing a better estimation of 
Doppler to be obtained. 

While some examples of the application of the invention have been described, it 
will be appreciated that the methods of the present invention can be applied 
widely to acoustic sounding and that many alterations and additions can be made 
without departing from the scope of the invention as defined by the following 
claims. 



